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Abstract 

We study the equation of state (EOS) of quark matter at zero temperature, 
using the Color Dielectric Model (CDM) to describe confinement. Sensible re- 
sults are obtained in the version of the CDM for which confinement is imposed 
smoothly. The two-phases version of the model turns out to give unrealistic 
results for the EOS. Chiral symmetry plays a marginal role and the quarks 
are massive till high densities. The deconfinement phase transition is smooth 
and unlikely to be first order. Instabilities of the quark matter and the gap 
equation are discussed. 

24.10.Jv Relativistic models, 

24.85.+p Quarks, gluons and QCD in nuclei and nuclear processes 



1 



Typeset using REVTeX 



I. INTRODUCTION 



The study of the Equation Of State (EOS) of Quark Matter (QM) has become a fash- 
ionable topic in view of the next experiments using heavy ions in program at RHIC(BNL) 
and at LHC(CERN) [|TJ. Furthermore the inner structure of neutron stars is now under in- 
vestigation: the connection between the composition of the star and the cooling time, which 
can be measured, allows to discriminate among the various models, indicating the possible 
existence of a quark matter phase (see, e.g. [0). The study of the equation of state of matter 
at high densities can also give usefull informations to traditional nuclear physics, since one 
can search heavy nuclei for precursor phenomena, both of the deconfmement and/or of the 
chiral restoration phase transition (for a review see [§]). 

In many model calculations of the deconfmement phase transition the frame of the MIT 
bag model has been used ||||. In such a way a first order deconfmement phase transition is 
obtained (apart from specific, ad hoc choices of the model parameters) and the deconfmement 
phase transition coincide with the chiral restauration one. At densities and temperatures 
slightly bigger than the critical ones the right degrees of freedom are already quarks, having 
current masses, and perturbative gluons. There are anyway indications, from lattice calcu- 
lations, that at temperatures bigger than the critical one non perturbative effects are still 
present in the quark-gluon plasma 0. 

In this paper we study the EOS of QM using the Color Dielectric Model (CDM) to 
describe confinement [0]. We shortly review CDM in sect. 2. This model has been widely 
used to study both the static and the dynamical properties of the nucleon. Morover it can 
be used to describe many-nucleon systems: for a two nucleon system it allows to compute 
a nucleon-nucleon potential qualitatively similar to the ones used in nuclear physics ||; in 
the case of a homogeneous, infinite system of nucleons the CDM can be used to construct a 
nonlinear version of the Walecka model ||. 

The aim of our work is to extend previous calculations of the deconfmement phase 
transition, where the same model has been used ||11|| . An important point of our calculation 
will be to fix the model parameters in order to reproduce the basic static properties of 
the single nucleon, as was already done in the study of the nucleon structure functions 
||12|| . We will later use the same parameters to study the EOS of QM, which we define 
as a system of totally deconfined quarks. In such a way the study of the QM's EOS will 
turn out to be a severe test for the different versions of the CDM, and we will be able to 
make some predictions of the properties of matter at high densities. Within CDM (with a 
double minimum potential for the scalar field), we will investigate the possibility of getting 
a scenario similar to the one described by the MIT bag model, with two phases undergoing 
a sharp first order phase transition (sect. 3). Our results show that such a description is 
incompatible with the CDM. We than study another version of the CDM (with a single 
minimum potential for the scalar field), where confinement is imposed more smoothly, and 
we get a sensible EOS for QM, without a sharp deconfmement transition. 

An important feature of this deconfined quark matter is that quark's masses are big till 
high densities. Chiral restauration and deconfmement do not occur at the same density. 
The reasons why chiral simmetry is restaured so slowly are discussed in sect. 4. 

In the last sections we analyze the properties of QM, as described in the CDM. In sect. 
5 the stability of quark matter is studied, using the technology of the response function. In 
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sect. 6 we consider a gap equation, trying to understand the formation of quarks' clusters. 
Finally sect. 7 is devoted to the concluding remarks. 



II. THE COLOR DIELECTRIC MODEL 



A. The model lagrangian 

In this section we shortly summarize the main features of the CDM. For a comprehensive 
review see the ones by Pirner |7| and by Birse [10]. 



We will use a chiral invariant version of the CDM, as the one used in Ref. []I1J and also 
in Ref. \T2\. The Lagrangian reads 



9 y 



C = ii/jj^daip H ip (a + ry 5 r ■ 7?) ip 

X 



-{d,xY -U(x) + ^{d,af + - U (a, vr) 



(1) 



where U(a, 7?) is the usual mexican-hat potential, as in ref. [fLS]] . £ describes a system of 
interacting quarks, pions, sigmas and a scalar-isoscalar chiral singlet field Xi whose potential 
U (x) has an absolute minimum for x = 0. In such a way, in the case of the single nucleon 
problem, the quarks' effective mass —gcr/x diverges outside the nucleon. The simplest 
potential for x is a quadratic one: 



U( X ) = \m\ 2 . 



(2) 



Using this potential, x fluctuates around zero, reaching this value asymptotically, at large 
distances. In the following we will refer to this version of the model as the Single Minimum 
(SM) one. Another possibility is to use a potential for x having a Double Minimum (DM): 



U( X ) = ^M\ 2 



1 + 



- 2 



X 
7M 




(3) 



where the absolute minimum is still in the origin and the relative one is in x — jM. Choosing 
appropriately the parameters it is possible to obtain solutions for the single nucleon where 
the x field interpolates between the relative minimum at the center of the nucleon and the 
absolute minimum at big distances. In the following, we will consider for DM only solutions 
where x h &s the above discussed behaviour. A complete analysis of the various possible 
solutions of the model 
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(in a non-chiral-invariant version) can be found in Ref. 
In Fig.|I] we show two typical solutions of the model for the single nucleon problem. 
As it appears, in the DM case the transition between the interior of the nucleon and the 
external region is sharper than in the SM, where all the fields have a very smooth behaviour. 
Correspondingly, the kinetic energy contribution (which comes mainly from the quarks) will 
be bigger in DM than in SM. This point will be relevant when studying the EOS of QM. 

The Lagrangian is chiral invariant, and it can be considered a confining version of the 
traditional a-model. An important point concern the value of the chiral fields, the pion 
and the sigma, in this model, in the single nucleon case. These fields are always near their 
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vacuum value. This point has been checked out numerically several times and some euristic 
explanations have been proposed [p~0, 14 . As a consequence, the chiral fields cannot 'wind' 
around the mexican hat potential. The pion is thus just a perturbation and cannot develop 
a non-trivial topology. 



B. Fixing the parameters 

The parameters of the model are: the chiral meson masses = 0.14 GeV, m a = 1.2 
GeV, the pion decay constant /„- = 0.093 GeV, the coupling constant g, and the parameters 
appearing in the quadratic (||) or in quartic \~ potential @. 

The free parameters are fixed to reproduce the basic properties of the nucleon. At the 
mean field level we use the hedgehog ansatz which is an eigenstate of the so called Grand 
Spin G = S + I , and is a superposition of various bare nucleon and delta states. 

In order to describe the single nucleon state we performed a double projection on linear 
and angular momentum eigenstates from the hedgehog, whose details can be found in Ref. 







For the DM version of the model we use the sets of parameters of Ref. |]T6| , for which a 
good description of the static properties of the nucleon was obtained. 

For the SM version such a set of parameters was not available in the literature. We 
fixed the parameters g and M to reproduce the experimental value of the average mass of 
the nucleon and of the delta, and the isoscalar radius of the nucleon. Choosing g = 0.02 
GeV and M = 1.7 GeV we got: (E N + E A )/2 = 1.112 GeV (exp.val.= 1.085GeV) and 
< r N >l isoscaiar = 0-82 fm (exp.val.= 0.79 fm). These values depend essentially only on the 
quantity G = \fg~M . A detailed presentation of the single-nucleon properties in SM version 



of the CDM will be presented elsewhere [17]. 



To perform an exhaustive analysis of the various versions of the CDM, we considered also 
the possibility of having an effective mass term for the quarks in which the \ field appears 
with a power different from one: m q = —g p a/x p - Studying this possibility in the SM version, 
and for p = 2, a good description of the single nucleon properties can be achieved using 
g = 0.02 GeV and M = 1.10 GeV ((E N + E A )/2 = 1.102GeV and < r 2 N >]£ acalar = 0.78fm). 

Several important differences exist between the CDM and the MIT bag model. In the 
latter model, inside the bag the quarks have current masses of few MeV. The bag is stabilized 
through the introduction of a big vacuum pressure, of the order of 150MeV/fm 3 . Perturbative 
gluons are considered to be the right degrees of freedom inside the bag. In the CDM model, 
in all versions, the effective quark mass is everywhere bigger than a number of the order 
of lOOMeV, hence chiral symmetry is broken and Goldstone bosons are the right degrees of 
freedom. We will come back later to this point, comparing the EOS of QM as computed in 
the chiral CDM with the one computed in a non-chiral version [18]. In the CDM a vacuum 



pressure is also present, coming from the x field: this pressure is roughly constant inside the 
nucleon in the DM case, where it equals U(x = jM) = M 4 ^ 4 , whilst in the SM case the 
pressure |M 2 x(r) 2 depends on r. It is important to stress that in the DM version of the 
CDM model the pressure is very small, of the order of few MeV, and this point will also be 
important when discussing the EOS of QM. 
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III. MEAN FIELD APPROXIMATION TO THE EOS OF QM 



In this section we will study the EOS of QM, using the Lagrangian of the chiral CDM 
(1-3). 

For QM we mean a system of totally deconfined quarks, described using plane waves, 
with the Xi pion and sigma fields having a constant value, given by the Euler-Lagrange 
equations. 

The total energy of QM in the mean field approximation is the following: 



E QM = 12V J ^ ] j^ + ( 9 jyd(k F -k) + VU(x) + VU(a,n = 0), (4) 

where kp is the Fermi momentum of quarks, x an d o are the solutions of the coupled 
equations 

dU(x) -Ps(x,°) /,n 

— 5 — _= -go- — — 2 — , (5) 

dx x=x x 
dU(a,n = 0) Ps(x,v) 

-= ~9 = , 6 

da °"=°" x 

and the scalar density ps(x, a) is given by 

p s (x, a) =< ^ >= 12 / ^ J(k F - k). (7) 



In the mean field approximation, for an homogeneous infinite system, < ff >= 0, so the 
pionic field is not contributing (of course it indirectly enters the EOS, because our model 
parameters are fixed in the single nucleon problem, where < 7r >^ 0). This is drawback 
of the mean field approximation, when applied to homogeneous infinite systems. A way to 



circumvent this problem is discussed by Ghosh and Phatak [[HJ. 

It is important to analyze the behaviour of the x field, both in the SM and in the 
DM version of the model, when the CDM is used to describe a collection of nucleons at 
increasing densities. First of all, in the DM model a critical density exist, for which the 
X field undergo a discontinous jump. In the appendix a proof of this statement is given, 
based on the study of a two-nucleon system for various internucleon distances. For distances 
smaller than a critical one, the x field, in the region between the two nucleon, will cease to 
interpolate between the two minima, as discussed in Sec.2B, and will stay near the relative 
minimum. This transition cannot be made continuous. In the DM version of the CDM, the 
deconfinement phase transition is therefore a discontinous first order transition. In the SM 
version of the model nothing similar can happen. Of course this is not enaugh to conclude 
that, in this case, the deconfinement phase transition is not first order. We will study more 
in detail what happens in Sec. 5 and 6. 

We will now compare the EOS of QM, as computed in our model, with the EOS of 
nuclear matter as obtained in the Walecka model PU|. In Fig.|2| our results for the EOS are 
shown. 



In the DM case we used the parameter sets 15,16 : 
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i): g = 0.059 GeV, M = 1.4 GeV,7 = 0.04,7? = 0.06 



ii) : g = 0.029 GeV, M = 1.2 GeV,7 = 0.06,7? = 0.06 

iii) : g = 0.0235 GeV, M = 1.6 GeV,7 = 0.03,7? = 0.06 

As it appears, the energy per baryon number in DM is very small, its value being below the 
one given by the Walecka model for almost all densities. The DM version of the CDM is 
therefore unrealistic when used to describe the EOS of QM. A similar result was obtained in 
a calculation performed in a non-chiral version of CDM [|ll|. The reasons for such a debacle 
are to be find in the very small value of the pressure in DM. If one would add 'by hand' a 
pressure's contribution of the order of lbOMeV/ p to the energy per baryon shown in Fig.|2| 
for DM, one will get a result similar to the one obtained using the MIT bag model. The 
problem is that such a big pressure cannot be obtained in the CDM, because corresponds 
to a solution for the single nucleon problem where the quark fields are very steep: as it has 
been shown by Leech and Birse lf21|j , the center of mass motion cannot be projected out 
consistently in this case, and the mean field approximation is no more a good starting point. 

In the SM with p = 1 case (see Sect.2B), on the other hand, the mean field approximation 
to the EOS of QM gives a sensible result: the energy per baryon number is bigger in the QM 
phase than in the hadronic phase for all densities smaller or of the order of p eq = 0.17 N/ fm 3 , 
the equilibrium density of nuclear matter. After this density the equation of state of QM 
and the one of nuclear matter seems almost equivalent, till densities of the order of 2p eq , 
after which the energy of QM is smaller than the energy of nuclear matter. We would like 
to remind that we have not modified the parameters of the model, but we are sticking to 
the ones fixed to the static properties of the nucleon, as discussed in Sec. 2. 

The last possibility we have considered is SM with p = 2. Also in this case, the EOS of 
quark matter that one obtains after fixing the parameters is too low, as it can be seen from 
Fig-0. 

We can conclude from analysis of all the versions of the CDM, that the most realistic 
EOS of the quark matter is obtained using the model in which confinement is imposed in 
the smoothest way, i.e. SM with p — 1. In the following sections, only this version will be 
considered. 

The most relevant feature of the result for SM p = 1, is the wide range of densities 
for which the EOS of traditional nuclear matter and the EOS of quark matter are almost 
equivalent. The difference in energy between the two phases is of some tens' MeV, only. It 
is remarkable that this 'almost equivalence' starts at a density of the order of p eq . A natural 
interpretation of this result is that in the case of heavy nuclei, some precursor phenomena 
of deconfinement could be seen (as swelling, for instance), but no dramatic change is going 
to happen in the system till much higher densities. In Sec. 6 we will analyze the possibility 
that the small energy gap between the two phases can be explained taking into account 
correlations in QM, using Bethe-Goldstone equation. 



IV. QUARKS' MASS AND CHIRAL SYMMETRY RESTAURATION 

We discuss now more in details the dependence of quarks' effective mass m q = —gcr/x 
on the density. In this model two different mechanisms are at work to reduce m q : one comes 
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from the chiral field sigma, that moves from its vacuum value —f w as the density increases 
(in the following we will not take into account a chiral-symmetry breaking term giving mass 
to the pion, because it is irrelevant to what we want to discuss). The other mechanism that 
modifies m q is confinement, through the \ field, which moves away from zero. As we will 
see, this second mechanism is the relevant one till high densities. 

We restrict our discussion to the SM p — 1 version of the model, the one which gives 
sensible results for the EOS of QM. In this case, eq.(5) can be formally solved, giving 

X = (-jjfs* < W> >) 1 / 3 . (8) 

This is not really the solution of eq.(5), because the scalar density ps =< 4>4> > still depends 
on x through the fermions' mass. It can nevertheless be used as an approximation to the 
real solution, if one neglects the difference between p =< tfj^ip > and p$ (this approximation 
is not too bad, because, as we shall see, quarks' masses are not decreasing very fast). 

The 'Maxican hat' potential, entering the Lagrangian (1), and parametrized as in Ref. 
is 
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U^) = ^ 2 + * 2 -fl?i (9) 

where m a is the mass of the sigma field. Using the formal solution for \ given in eq.(8), one 
can compute the shift A of the sigma field from its vacuum value 

a=-U + A (10) 



1/3 

(^) 2 / 3 EE C(^) 2/3 . (11) 

Substituting in eq.(ll) typical numbers, and in particular the ones we used in Sec. 3, the 
dimensional coefficient C in eq.(ll) turns out to be small, of the order of 0.1 fm, and A is 
therefore also small, if compared to /„., till very high densities. 

In Fig.|3| we compare the effective quark mass, as computed taking into account both 
the reductions coming from the \ field and from the a field (solid line), with the effective 
mass when the a field is kept fixed at — f n (dashed line). At a density of order of 3p eq , the 
effective mass of the quarks is reduced from its value at p eq by a factor ~ 0.8. Comparing 
the two lines of Fig.||, one can see that most of the effect is due to the confining \ field. 



A 
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V. QUARK MATTER INSTABILITIES 



The more direct way to study the instabilities of a many-body system is to look where its 
compressibility K becomes infinite. The compressibility is related to the pressure through 
the following relations 

P = fW10 (12) 
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(13) 



Here (E/N) is the energy per particle and p is the density. 

In Fig}4| the compression modulus K~ l is shown (here and in the following we will refer 
to the SM version of the CDM, only). As it can be seen, quark matter becomes unstable, 
in our calculation, at a density smaller than pi nst ~ p e g/2. This instability is not due to the 
fact that in this region the energy per particle of quark matter is higher than that of nuclear 
matter, but it is related to the possibility of creating undamped density fluctuations in the 
system, without spending energy. 

To clarify even more this point, one can try to reproduce the same result, studying 
the collective excitations of the system, which can be found searching for the poles of the 
propagators of the scalar fields, dressed by the particle-hole polarization propagator. In the 
following we will only discuss the collective states coming from the propagation of the x an d 
the sigma field, leaving the discussion of the pion propagator, and the related phenomena, 
to a future study. 

To define the propagators of the x an d the sigma, we expand these fields around their 
mean field value: 

a = a + a (14) 



X = X + X- (15) 

The mass term reads thus (< n >= 0) 

Tpaip ip(a + a)if) ipaip ipaip ipaip 
9 =9 ( - , ~x = 9^— + 9— 9—^-X + - (16) 

x ix + x) X X X 

where we have expanded the x field till first order in the fluctuation. Since we are looking for 
the instabilities of the system, i.e. for the situations in which an (arbitrary small) fluctuation 
around the mean field develops spontaneously and propagates undamped, it is enough to 
consider a first order expansion in \. 

From the linearized lagrangian (|TJj), we read the couplings between the scalar fields x 
and a and the quark fields 

9x = -9^ (17) 



9 

9a = -• 

X 



(18) 



The masses corresponding to the x an d to the a fluctuations are given by the following 
relations 



<9 2 

M* 2 =< — 
ox 



u{x)-g 



X 



>= M 2 - 2g^$- = 3M 2 



(19) 
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(20) 



In the previous equations, the brakets <> indicate mean value in the mean field approxima- 
tion, where the fluctuations are set equal to zero. In eq.(^), we have used the field equation 
for x @ ; to eliminate the dependence on the fields. The mass of the a is only slightly reduced 
from its value at zero-density, because of the slowness of the chiral symmetry restauration 
(see the discussion in the previous section). 

The mean field propagators of the scalar fields are thus 

« = i/( g ;- m f + ^) (2i) 



We define the scalar polarization propagator in the following way 

d A k 



n.(?) 



(27T) 



rTr[G(k)G(k + q)], 



(22) 



(23) 



where G(k) is the single quark relativistic propagator in the medium which can be decom- 
posed in the Feynman propagator and a density dependent correction 



G(k) = ( 7 % + m q ) 



1 



VK 



m q + irj 



E*(k) 



5{k Q - E*(k))9{k F - \k\) 



= G F {k) + G D {k). 

(24) 



The effective quark mass is computed from mean field equations: m q = —gu/x-, and E*(k) 



m 2 + k 2 . For a review of the formalism, see Reff. [ P0| , P2|] . 
In accordance with the mean field approximation, we neglect the vacuum fluctuations 
effects, and we take only the density-dependent part of the polarization propagator |2"2fl . In 
such a way we avoid possible complications coming from the ill-defined zero density vacuum 
state. An explicit expression for Il s can be found in Ref. p2" |. 

We now consider the propagators of the scalar fields, modified by the polarization propa- 
gator insertion, i.e. the propagators in the RPA approximation. The x an d the a propagators 
mix up in the RPA approximation. They are the solution of the following coupled equations 

D™ = Dq + Dq g 9 Il s (g a D™ + g^D**) 
D xx = D x + D x ^ U s (g^D xx + g,D° x ) 
D^ = D x g^Yl s (g d D™ + g^) 

D ax = Dq ga n s ( g % D™ + g B D°*). (25) 

In the previous equation, the dressed propagators correspond to various situations in which 
the two mean field propagators Dq and D x mix among themselves via the polarization 
propagator. These propagators can be decoupled, giving the following RPA propagators 
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D^ = D^l-g^U s g^)/e 
D™ = D x (l-g a Tl s g a D°)/e 
DX° = IJ$gzIL a g a DZ/e 
D^ = D^g a U sgjc D oX /e 

e = 1 - ( 9ji U s 9jc D x + g a U s g a D%). (26) 

The collective excitations of the system correspond to the zeros of the denominator e(g M ). 
Instabilities are the collective states at zero energy transferred [22], i.e. the solutions of the 
equation 

e(q = 0,q) = (27) 

The solutions of this equation are shown in Fig.[|, as a function of the transferred momentum 
and of the density of the system. In the region enclosed by the line, corresponding to 
the solutions of eq. Q2"7|), the system is instable; in particular, when these instabilities 
exist for zero transferred momentum the system develops spontaneously undamped density 
fluctuations. The range of densities for which the system is found to be instable (at zero 
momentum transferred) coincides with the range found studying the compressibility. 

Concerning the interplay between chiral symmetry and confinement in the development 
of the instabilities, we found again that chiral symmetry plays a marginal role in the model: 
if one takes into account the x propagator only, the instability region shown in Fig.|4] is 
almost unchanged. 



VI. TWO BODY CORRELATIONS: GAP EQUATION 



We will now study two-body correlations, using the formalism of the gap equation p3 
The gap is the difference between the energy of the incorrelated pair and the energy of the 
correlated pair. It is thus positive if the potential is attractive. We will use as the residual 
interaction the one arising from the exchange of a x or a a. Of course we are able to take 
into account only two-body correlations, and three body-correlations can be even more 
important. 

The gap equation reads P~~ 



l£< k ,- kW -W >M ^ (28) 



Here we are following the convenction [23] that an attractive potential is positive. The 



potential V is evaluated between incorrelated states, described by plane waves. is the 



single-particle energy, mesured relative to the chemical potential /i = ^Jk 2 F + m 2 . We will 
not use the relativistic reduction of because the relativistic corrections are not totally 



negligible. Therefore = ^jk 2 + m 2 — \jk 2 F + m 2 , where m q is the quark mass as computed 
in the mean field approximation. 

As we have already said, the potential V arises from the exchange of the scalar fields' 
fluctuations x an d cr. The masses of this fluctuations and the couplings between the fluc- 
tuations and the quarks' fields, have been obtained in the previous section, expanding the 
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lagrangian around the mean field approximation. Of course this expansion, which was suffi- 
cient to study the instabilities of the system, is in general quite a crude approximation. The 
results obtained for the gap are thus only the first step in the study of the difficult problem 
of the clusterization. 

The potential V is a Yukawa potential, but an extra factor m 2 q / E* (k)E* (k 1 ) has to be 
included (E* (k) = ^Jk 2 + m 2 ), because of the choosen normalization of the quarks' spinors. 

Following Ref. since the resulting gap A k is much smaller than the Fermi energy, 
the integrand appearing in eq.(p8|) is sharply peaked near £ = 0, and then A k ~ A kp = A, 
i.e. the gap is almost independent on the momentum. 

We show in Fig.(| the resulting gap. It shows a strong dependence on the density, and 
for densities of the order of p eq is already totally negligible. 

The gap equation can be solved analytically for small value of the gap. The result is 

A « 8^exp( * k l{*™ q ), (29) 

2m q yy <k F \V\fa F > h 1 ; 

where the matrix element of the potential is 

< k F \V\4>k F >= k F I sinkpxV (x)smk F xdx. (30) 
Jo 

Since our residual interaction V is always attractive, it doesn't exist a critical density at 
which the gap is exactly zero, and this value is reached only asymptotically. 



VII. CONCLUSIONS 

In this paper we have studied quark matter, using the non-perturbative tool of the CDM 
model. Let us summarize our main results: 

- of all the considered versions of the model, only one gives sensible results, i.e. the one 
in which confinement is imposed in the smoothest way. The other versions of the model give 
an exceedingly low energy per baryon number for the quark matter. 

- the SM (p=l) version gives an EOS for the quark matter which is almost identical to 
the EOS of nuclear matter as computed using the Walecka model, for the range of densities 
Peq < P < 2 p eq . In this range, the difference in the energy per baryon number between 
the nuclear matter and the quark matter is very small. Taking into account the theoretical 
incertitude in the fixing of the parameters (Sec.2B), this energy difference is of the order 
of 20 MeV. At densities smaller than p eq the energy difference rapidly increases, and for 
densities higher than 2p eq the quark matter is the energetically most favourable state. An 
important point is that the minimum for the EOS of quark matter is at a density of the 
order of p eq , and this result does not depend on the fine tuning of the parameters. 

- the mass of the quarks remains big (of the order of 100 MeV) till high densities, much 
higher than the density at which quark matter becomes the ground state. The deconfinement 
phase transition and the chiral symmetry restauration arise at totally different densities. 

- the quark matter (in SM and p=l) becomes unstable at low densities, of the order of 
p e g/2. The instability can be obtained both from the study of the compressibility (where the 
compressibility becomes negative the system is unstable) and from the study of the collective 
sti cites ciij zero energy transfer. The two method give the same critical density. 
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- the process of clusterization can be studied considering the correlations between the 
particles, beyond the mean field approximation. The gap we obtained seems rather small. 
We have to bear in mind that we have oversimplified the problem, by linearizing the resid- 
ual interaction (and thus getting an approximate propagator for x, good only for small 
fluctuations), and by considering only two body correlations, where the three body ones are 
probably the most relevant. 

To conclude, we would like to consider three possible applications of the model. 

- Cooling of neutron stars (see C.J.Pethick in 0). 

A mechanism called URCA has been invoked to explain the rapid cooling of neutron stars. 
This mechanism proceeds via the exchange of electrons between neutrons and protons, which 
cool down emitting neutrinos and antineutrinos: 

n — > p + e~ + u e (31) 

p + e~ — > n + u e (32) 

A minimal fraction of protons is required, in order to fulfil momentum and energy conser- 
vation. This critical fraction is of the order of 1/9. Using traditional nuclear physics models 
to compute the protons' fraction, one gets numbers slightly smaller than the critical one, 
and the URCA mechanism cannot start. 

Another possibility is to invoke the presence of quark matter in the core of the star, 
and to consider reactions similar to the one previously described, but with the electron now 
exchanged between up and down quarks. In this case the problem is that, considering quark 
matter as described by the MIT bag model, quarks are massless and the phase space is 
thus zero. Therefore one would need a massive quark matter phase, and the possibility of 
reaching this phase at the density of the core of neutron stars, tipically of the order of 5p eq . 
This situation is actually the one described by the SM (p=l) version of the CDM. URCA 
mechanism should therefore be possible, and with an high luminosity, too. 

- Energy released in supernova explosion. 

Using a traditional nuclear physics approach, the energy released in supernova explosion 
is generally too small. A softer EOS could solve the problem, but if one uses e.g. the 
MIT bag to study matter at high density, the deconfinement phase transition is reached 
at densities larger than the one presumably reached in the collapse of the star. The EOS 
for matter at high density as computed in the CDM, is softer than the EOS of nuclear 
matter, and presumably a similar result will be obtained when computing neutron matter. 
Furthermore the softening starts at densities of the order of 2p eq . 

- EMC effect and swelling of the nucleon. 

To conclude let us consider the problem of the possible swelling of the nucleons embedded 
in a nucleus. If one considers,e.g. electron-scattering on heavy nuclei, one realizes that the 
swelling is a sensible mechanism, but it must be of the order of some 5% in order to be 
realistic. The real problem is thus not the one to obtain a swelling, but to obtain a not too 
big effect. In other words, the nucleons have not to dissolve when embedded in a nucleus. 
Since the minimum of the EOS of the quark matter is for a density near p eq , only in the 
center of heavy nuclei some swelling mechanism can appear. The exact amount of swelling 
depends on the precise difference in energy between the quark matter and the nuclear matter 
at densities ~ p eq , and is beyond the possibility of the present calculation. 

We are now carrying out researches in all the directions previuosly outlined. 



12 



APPENDIX: 



In this appendix we discuss the behaviour of the \ field in the DM version of the CDM, 
for a two nucleon system, as a function of their distance d, which is related to the density 
of the system p oc d~ 3 . We will show that a critical internucleon distance exists, at which 
the x field has a discontinuous behaviour. 

For simplicity we will omit the chiral fields, the pion and the sigma, from our discussion, 
considering the following Lagrangian: 



where the potential for x is the one given in eq.(3) and shown in Fig. 5a. We have labelled by 
Xm the relative minimum and by Xi the other value of x f° r which U(xi) = U(x m ) = M 4 rj 4 . 

We consider a system made of two clusters of three quarks each, with an intercluster 
distance d. In the mean field approximation the total energy of the two-nucleon system is 
given by 



where e q is the energy of the single quark, T x is the energy of the x field associated with its 
spatial fluctuations Q and E x is the energy coming from the potential U(x)- 

When the intercluster distance d is large, the x field interpolates between a number 
slightly bigger than Xm and a small number, smaller than xi, i n the internucleon space (see 
Fig. 5b, the solid and dotted line). 

When the distance d is reduced, the value of the scalar field in the internucleon region 
increases. We compare now two possible solutions: one in which the x field has a minimum 
value in the internucleon region, equal to Xi> an d a second in which the x field remains 
almost constant in the internucleon region (solid and dashed line in Fig.5b). The second 
solution has a smaller energy than the first one. In fact: e q is smaller, because the quarks 
move in a single big well, instead of moving in a double well; T x is smaller because x is 
fluctuating less; E x is smaller because x is n °t moving (in the internucleon region) through 
the relative maximum of U(x)- Since the second solution has a smaller energy, the value 
of x m the internucleon region will not smoothly increse, as the distance d is decreased, 
but will jump, at a certain critical distance, from the behaviour described by the solid and 
dotted line in Fig. 5b to the one described by the solid and dashed line. The internucleon 
value of x can be assumed as an order parameter that undergo a discontinuous change as 
the density is increased, thus the transition is first order. 



This energy is not the kinetic energy of the x field, which is zero, because x is a scalar field and 
it is assumed to be time independent. 



c = i^d^ - §m + -{d,xf -u(x), 

2 



(Al) 



E 2 = 6e q + T X + E X 



(A2) 
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FIG. 1. Typical solutions of the model for single nucleon problem in the SM (a,b) and DM 
(c,d) versions. 
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FIG. 2. EOS of QM and EOS of nuclear matter in the Walecka model (dotted line). We show 
the EOS of QM in CDM for the DM version of the model (dashed lines i) ii) hi) ) and for the SM 
version with p=l (solid line) and p=2 (dot-dashed line) . 
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